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We present energy filtered electron emission spectromicroscopy with spatial and wave- vector reso- 
lution on few layer epitaxial graphene on SiC(OOOl) grown by furnace annealing. Low energy electron 
microscopy shows that more than 80% of the sample is covered by 2-3 graphene layers. Cls spec- 
tromicroscopy provides an independent measurement of the graphene thickness distribution map. 
The work function, measured by photoelectron emission microscopy (PEEM), varies across the sur- 
face from 4.34 to 4.50eV according to both the graphene thickness and the graphene-SiC interface 
chemical state. At least two SiC surface chemical states (i.e., two different SiC surface structures) 
are present at the graphene/SiC interface. Charge transfer occurs at each graphene/SiC interface. 
K-space PEEM gives 3D maps of the |feii | n — n* band dispersion in micron scale regions show that 
the Dirac point shifts as a function of graphene thickness. Novel Bragg diffraction of the Dirac 
cones via the superlattice formed by the commensurately rotated graphene sheets is observed. The 
experiments underline the importance of lateral and spectroscopic resolution on the scale of future 
electronic devices in order to precisely characterize the transport properties and band alignments. 

PACS numbers: 73.22.Pr, 61.48.Gh, 79.60.-i 
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I. INTRODUCTION 



With the demonstration of GHz FETs based on epi- 
taxial graphene grown on SiCji^ this material has be- 
come the leading candidate for graphene-based electron- 
ics. Nonetheless, exploiting the remarkable properties of 
graphene for carbon based electronics remains an impor- 
tant challenge. The band structure and transport prop- 
erties of graphene must be either preserved or modified 
in a reproducible fashion on typical device scales. While 
a good deal of work has already focused on monolayer 
graphene grown on the SiC(OOOl) (Si-face),— the ability 
to grow thin graphene films on the SiC (0001) (C-face) 
has only recently been demonstrated^ C-face films offer 
a particular advantage because of their rotational stack- 
ing that effectively decouples adjacent graphene layers^ 



This leads to very high mobilities^- and allows devices to 
be less sensitive to thickness variations. While these sys- 
tems continue to make progress towards realistic carbon 
electronics, significant research problems remain. One in 
particular is the study and control of the graphene-SiC 
interface. 

When graphene is grown on the Si-terminated SiC 
on SiC(OOOl) (Si-face), the first graphene layer grows 
on an insulating carbon buffer layer with a (6-\/3 x 
6\/3)R30 o symmetry. 9 This buffer layer has a graphene 
structure and can be isolated from the SiC by inter- 
calating hydrogen between the SiC substrate and the 
buffer layer ji^ However, mobilities in this isolated buffer 
layer remain low suggesting either some type of prior 
disorder in the layer before hydrogenation or an effect 
caused by the hydrogenation itself. Very little is known 
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about the C-terminated SiC(OOOl) graphene-SiC inter- 
face. However, it is known that the C-face and Si-face 
interfaces must be very differentjii X-ray studies show 
that the atomic density gradient at the interface is differ- 
ent for Si-face and C-face) 12 ' 13 and core level photoelec- 
tron spectroscopy shows clear differences between the two 
interfaces^ Also, unlike the Si-face, the C-face interface 
is known to have two coexisting structures (at least in 
the early growth phase). Scanning tunneling microscopy 
(STM) studies have shown that, in UHV growth condi- 
tions, poorly ordered (2 x 2) and (3 x 3) surface recon- 
structions exist below the first graphene layerJ^. How- 
ever, these structures may disorder or simply not exist 
when growth occurs at higher temperatures. In high tem- 
perature furnace growth, there is no real evidence that 
the interface has a reconstruction. It is either an ordered 
(lxl) or a disordered reconstruction. The most impor- 
tant observation demonstrating that the C-face interface 
is very different from the Si-face is that C-face graphene 
has a rotational stacking very different from that of Si- 
face graphene even when growth temperatures are the 
same^iS This implies that the registry forces, and thus 
the interfaces, must be different on the two surfaces. Fi- 
nally the rapid growth of C-face graphene at tempera- 
tures lower than those observed on the Si-face point to a 
significant difference in the chemistry of the interfaced 

In this work we present a detailed study of the C-face 
graphene-SiC interface. To date, most of the experimen- 
tal techniques used to study the atomic and electronic 
structure are area-averaged and are therefore not sensi- 
tive to variations on the micron scale. Near field methods 
such as scanning tunneling microscopy (STM) can reduce 
the probed area by two orders of magnitude. However 
STM does not probe the interface itself and still provides 
very little information on the intervening length scales, 
which arc precisely those of interest in many potential 
device applications. Recently this intermediate length 
scale has begun to be explored using Low energy elec- 
tron microscopy (LEEM). Luxmi et alJ£ have studied 
the morphology of both UHV and argon furnace growth 
C-face graphene. Their work revealed a great deal of 
graphene thickness variation in these thick argon grown 
films. 

We focus on the spatial variation of both the electronic 
structure and chemical bonding of C-face graphene-SiC 
interface. The studies were carried out on thin C-face 
graphene films grown by a controlled Si sublimation tech- 
nique. In particular, we investigate the chemical homo- 
geneity of the interface, and correlate it with changes in 
doping, the graphene work function and graphene's 2D 
band structure near the Fermi level. We show that the 
interface is very complicated with local chemical changes 
that are not all associated with the local graphene thick- 
ness. To carry out these studies, we use LEEM, photoe- 
mission electron microscopy (PEEM) , and X-ray photoc- 
mission electron microscopy (X-PEEM). In addition, by 
using a suitable lens configuration the focal (or diffrac- 
tion) plane of the PEEM can be imaged to give parallel 



momentum resolved d ispersion curves E^kx^ky ). This 
technique is known as fc-resolved photoemission electron 
microscopy (fc-PEEM). Imaging the focal plane in PEEM 
produces a map for all azimuths simultaneously. Conser- 
vation of the component of the electron wave vector par- 
allel to the sample surface automatically transforms this 
map into one of photoelectron intensity as a function of 
(k x , fc y ), that is a horizontal cut in reciprocal space. Com- 
bined with energy analysis, this produces an image in re- 
ciprocal space of the local intensity as a function of wave 
vector parallel to the surface. For example, the Fermi 
surface can be acquired in a single shot experiment^ 



II. EXPERIMENT 

The substrate used in these studies was a 6H conduct- 
ing SiC(OOOl) from Cree Inc. Before graphene growth, 
the sample was first H2 etched for 30min at 1400 °C . 
The sample was then grown in an enclosed graphite 
RF furnace using the confinement controlled sublima- 
tion process, CCS. 4 The growth was done at 1475 °C for 
20mins. 

Before all measurements, the sample was annealed at 
500 °C for 1 minute in UHV to remove surface contam- 
ination. The surface cleanliness was checked by Auger 
electron spectroscopy and X-ray photoemission (XPS). 
LEEM was used to quantify the graphene layer thick- 
ness. LEEM experiments were carried out using a com- 
mercial Elmitec PEEM/LEEM III with base pressure of 
4 x 10~ 8 Pa. The energy filtered X-PEEM experiments 
were conducted on TEMPO beam line of the SOLEIL 
synchrotron (Saint Aubin, France) using a NanoESCA 
X-PEEM (Omicron Nanotechnology GmbH), 1 ^^ A dou- 
ble pass hemispherical energy analyzer was used to com- 
pensate single analyzer aberrations. This resulted in a 
PEEM energy resolution of 0.2eV with a lateral resolu- 
tion of ~ lOOnm for core level emission. Experiments 
were conducted in Ultra High Vacuum (UHV) at a pres- 
sure of 6 x 10 -9 Pa. 

The X-PEEM image series was acquired over the 
photoemission threshold region and the Cls core level 
(/ii/=654.3eV). A 53/im filed of view (FoV) was used 
with 12 kV extraction voltage. For the real space PEEM 
mode (threshold and Cls), a contrast aperture of 70/im 
was used, the lateral resolution was estimated to be 100 
ran. The analyzer entrance slit was set to 0.5 mm with 
a lOOeV pass energy to give a resolution of 200 meV. 
The resolving power of the TEMPO beamline is approx- 
imately 5000, giving an overall estimated energy reso- 
lution better than 250 meV. Dark images were acquired 
with the MCP turned off in order to remove camera noise. 
Flat Field images were acquired to correct for MCP de- 
fects. The parabolic non-isochromaticity of the instru- 
ment was corrected for all images 

The /c-PEEM results were acquired using the same in- 
cident x-ray spot position on the sample as the X-PEEM 
analysis (beam size ~ 100/imx50/im provides uniform 




FIG. 1. LEED patterns from a nominally 3-layer C-face 
graphene sample with primary electron energy of (a) 76eV 
and (b) 126eV. 

illumination). Because of the high extraction voltage be- 
tween the sample and the lens, the wave vector resolution 
and the dimensions of the reciprocal space image are in- 
dependent of the photon energy for the typical spectral 
ranges used in these experiments. The lateral spatial res- 
olution in the fc-PEEM mode was purposefully reduced 
by operating with a full open aperture. This was required 
to image a sufficient portion of reciprocal space in order 
to cover a full Brillouin zone. In this setup, the area of in- 
terest on the sample surface is chosen by a field aperture 
situated in an intermediate image plane that was closed 
down to about 7/im. A transfer lens then projected the 
1500/im diameter disk of the focal plane via the energy 
analyzer onto the detector giving a 2D fc-space dimension 
of about ±2.5A around the T-point. The spectrometer 
resolution was 200 meV, the photon bandwidth 20 meV, 
and the wave vector resolution ^0.05A . The detector 
response was corrected by the flat field of the detector 
and camera defects were eliminated using dark images. 
Operating this way gave a ~ l/i lateral resolution for the 
band structure imaging. 

III. RESULTS 

While Si-face graphene is known to grow oriented 30° 
relative to the principle SiC(2130) direction, multilayer 
C-face graphene is known to be stacked with sheets 
within the stack having multiple rotation angles peaked 
at 30° and 0± ~ 7°. These C-face graphene rotations 
are due to interleaved rotated graphene sheets with non- 
Bernal stacking (i.e non-60° rotations). In thin C-face 
samples there are only a few rotations. This is demon- 
strated in the LEED patterns in Fig.ffl 

Because the films are thin, the LEED pattern clearly 
shows the six- fold SiC substrate spots (Si and 5*2 ). In 
addition, the LEED shows a set of graphene spots ro- 
tated by 30° with respect to the S2 (31 and #2) and 
a second set of spots rotated ^4 — 5° relative to S*2. 
Weaker graphene arcs around 0° are also visible. The 
SiC substrate LEED pattern is 1 x 1 and no additional 
diffraction spots that would suggest a significant recon- 
struction are visible. However, these LEED images are 
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FIG. 2. (a)-(c) Typical LEEM images (10/zm FoV) of the 
graphene surface as a function of electron energy (a) 3.0eV (b) 
4.8eV (c) Reflectivity spectra extracted along the dotted line 
in (b) showing clearly the variation in the number of intensity 
minima and therefore the number of graphene layers. 



macroscopic area-averaged results, and are not expected 
to be sensitive to microscopic variations in the interface 
structure. More detailed spatial information is obtained 
using LEEM and X-PEEM as discussed below. 



A. LEEM 

LEEM data were obtained by using the (0,0) specu- 
lar back-scattered electron beam. Figure [2] shows typical 
bright-field images with a FoV of 10/im for electron ener- 
gies of 3.0 and 4.8eV. A full image series was acquired by 
varying E from 1.5 to 12.9eV (using a 50meV step size). 
The low energy onset of the back-scattered electron signal 
depends on the potential just above the surface and thus 
on the local work function. Clear differences in contrast 
as a function of E are observed across the field of view 
suggesting distinct work function for different regions of 
the surface. 

In two-dimensional layered systems there are os- 
cillations in the LEEM reflectivity at low electron 
energiesj22r— Several groups have used these reflectiv- 
ity oscillations to determine the number of graphene lay- 
ers in the epitaxial film. Figure [2fc) shows the spectra 
along the dotted line in Fig. [2Kb). We observe oscilla- 
tions between 1.5 and about 7.5eV. In graphite there are 
band gaps below OeV and above 7eV along the TA direc- 
tion. In an ideal multilayer graphene film, the number 
of minima in the reflectivity between successive Bragg 
peaks gives directly the number of graphene layers. We 
have similarly extracted pixel-by-pixcl reflectivity curves 
from the image stack and mapped the LEEM reflectivity 
across the 10/im FoV and used the number of minima 
in each reflectivity curve to produce a thickness map of 
the C-face graphene film. Figure EKa) shows the typical 
reflectivity curves for each distinct contrast region in the 
FoV. Two slightly different curves without a clear oscil- 
lation between 1.5 and 7eV are observed. They are both 
attributed to the zeroth layer or C-terminated layer of the 
SiC substrate. This layer is thought to interact strongly 
with the SiC substrate through sp 3 bonding. Luxmi et 
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FIG. 3. (a) Typical reflectivity curves extracted from the 
LEEM dataset showing to 3 oscillations below 7.5eV. (b) 
Graphene thickness map (FoV=10/im) generated by counting 
the number of minima in the reflectivity curves. 



al. 16 have shown even flatter reflectivity curves for the 
zeroth layer (OML) for surfaces prepared at higher tem- 
perature under an argon back pressure. The thickness 
map constructed from the intensity minima is shown in 
Fig. |3jb) using the same color coding. The data shows 
that more than 80% of the surface is covered by 2- or 
3-layer graphene. 

The reflectivity oscillation can be understood as quan- 
tum interference between electrons reflected by differ- 
ent graphene layers^ Hibino et al£& pointed out that 
although the conduction band is continuous in bulk 
graphite between 4.3 and lleV along FA, few layer 
graphene should have discrete states, thus the reflectivity 
oscillations are correlated with the electron structure of 
the thin films. The quantized conduction band states en- 
hances the transmission of the incident electrons produc- 
ing the dips in the reflectivity curves. This is confirmed 
by the good agreement between the experimental minima 
positions and the resonant energies predicted by a tight- 
binding calculation. For an m-layer thin film the bulk 



> 

>- 
E? 

<D 
C 
UJ 



. □ Experiment 




• 


• Theory ^ 


f 






• 


" I 

E 




• 




1 2 


3 


4 



Graphene thickness (ML) 

FIG. 4. The position (boxes) of dips in reflectivity curves 
from Fig. [2fc) and Fig. [3ja), compared with theoretical tight 
binding estimation based on tight- binding model (•). The 
error bar on the experimental data indicates the correlations 
between the spread in the energy of the reflectivity minima. 



band dispersion has discrete energy states when the wave 
vector satisfies the quantized condition E = e— 2t cos Ka, 
where e is the energy of the band center, t is the transfer 
integral, and a is the interlayer distance. At these values, 
the dips in the reflectivity are predicted by projecting the 
wave vector onto the calculated band structure along T— A 
direction, the normal to the graphene layers. 

Figure 0] shows the comparison of experimental reflec- 
tivity minima with the discrete energy levels predicted by 
the tight-binding theory. The spread in the experimen- 
tal energy of the reflectivity minima in different regions 
of the FoV is indicated as an error bar. The centre en- 
ergy of the tight-binding theory calculation is at 3.6eV 
compared to 3eV used by Hibino for Si-face graphene. 20 
Note that the third calculated minimum (highest energy) 
in the nominally 3ML region is at a slightly higher en- 
ergy whereas the first and second minima are at slightly 
lower energy than the experimental data. This trend is 
also seen in the 2ML minima. The difference may be 
due to the approximations inherent in the tight binding 
calculation. 



B. X-PEEM 

While LEEM measurements give information on the 
structural spatial variation of the graphene thickness, 
lateral variations in the electronic structure of epitax- 
ial graphene have never been mapped. XPEEM offers a 
unique method to begin to understand the role the SiC- 
graphene interface plays in graphenes electronic struc- 
ture. In this section we present the first sub-micron 
chemical and electronic mapping of the graphene-SiC in- 
terface. 
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Work Function (eV) 

FIG. 5. (a) Threshold image at E — iJp=4.45eV with a 53/im field of view showing clear intensity contrast due to work function 
variation across the surface, (b) Work function map of the FoV obtained from a pixel by pixel analysis of the threshold spectra, 
(c) Histogram of the work function values in (b) showing five distinct Gaussian distributions, corresponding to five distinct 
work function values. 



1. Work function 

Probing the transition from the mirror reflection of 
the electrons to the back scattering regime, commonly 
referred to as mirror electron microscopy to low energy 
electron microscopy (MEM-LEEM) transition is highly 
sensitive to the local variations of the electric potential 
just above the surface, as small differences in the lat- 
ter determine large differences in the electron reflectiv- 
ity. In energy-filtered PEEM, at high photon energies 
(654. 3eV in these experiments), the photo-emitted inten- 
sity at threshold is directly related to the work function. 
An example of a raw image obtained in the threshold 
region is shown in Fig. [Sfa). As E — Ep is scanned, 
dark areas become bright and vice-versa, giving rise to 
spatial contrast as a function of photoelectron kinetic en- 
ergy. This is direct evidence for a distribution of work 
function values over the sample surface. 

After correction for the Schottky effect due to the high 
extractor field, Ai? = 98 meV for 12 kV, 23 the photoemis- 
sion threshold spectra can be used to directly measure the 
local work function. The threshold spectrum is extracted 
from each pixel in the FoV (pixel area 65 x 65nm 2 ). The 
position of the threshold is obtained from a fit using a 
complementary error function, 



1(E) = A • erfc[^-=— ] + 7 min . (1) 
v 2c 

where <I> is the work function and a the half- width of the 
rising side of the secondary electron peak (O.leV here). 
We note that because the theoretical shape of the photoe- 
mission onset is modeled, this method is a more reliable 
way of obtaining absolute work function values than sim- 
ply extrapolating a straight line down to zero intensity. 
The results of this analysis are presented in the form 
of a map of the work function within the field of view 
in Fig. EKb). A histogram of the work function values 



across the whole FoV is shown in Figure [5jc) . Five Gaus- 
sians, with width ±25meV, are able to describe the work 
function frequency distribution suggesting that there are 
only five distinct values of the work function. Actually, 
as we show when we discuss the XPEEM results, the 
peak in the work function distribution centered at 4.46eV 
includes two slightly different different work functions 
(A$ = lOmeV). While these two regions of the surface 
have nearly the same work function they can be distin- 
guished by their very different core level spectra. 

The work function distribution spans a range from 
4.34eV to 4.50eV; a range well below the work function 
of bulk graphite (4.6eV); thus our results confirm that no 
significant part of the surface within the FoV consists of 
many layer graphene or graphite. On the contrary, this 
epitaxial film is indeed near-uniform few layer graphene, 
with local variations in the graphene coverage. However, 
we cannot necessarily attribute each peak in the work 
function distribution to a distinct graphene thickness. 
For example, recent Kelvin force microscopy showed that 
the work function difference between one and two layer 
graphene is 0.135eV<21 The graphene thickness calibra- 
tion, based on the local Cls XPEEM spectra presented 
below, excludes such a difference between regions with 1 
and 2 layer graphene, therefore the observed variations 
cannot be due to only thickness changes. This is partic- 
ularly true for regions of the surface with work function 
values of 4.46-4.49eV. For these regions the Cls spec- 
tra (as discussed in the next section) indicate that the 
graphene film is in fact very thin. 

There are also more substantial changes to the thresh- 
old spectra than the work function measurements sug- 
gest. To demonstrate this, Fig. [Ha) shows two examples 
of threshold spectra for regions of the sample with a small 
work function difference of ~ 30meV (4.42 and 4.45eV). 
These regions will be analyzed in more detail using both 
core level XPEEM and the fc-PEEM data in the next 
section. There are interesting variations in the structure 
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FIG. 6. (a) Local threshold spectra extracted from regions in 
Fig. Ob) with work functions of 4.42 and 4.45eV. The best 
complementary error function fits of the rising edge of the 
photoemission threshold are indicated by the solid line, (b) 
Second derivative smoothed threshold spectra from (a) show- 
ing peaks corresponding to the empty conduction band struc- 
ture. 



of the secondary electron (SE) peak up to 4-5eV above 
the vacuum level (i.e. above the photoemission thresh- 
old). The SE peak structure can be related to the con- 
duction bandit Figure [f^b) shows the negative of the 
second derivative threshold spectra. The data have been 
smoothed in order to more clearly see the peak structure. 
The low work function regions in Fig. [5jc) have a single 
main structure in the SE at ~ 6.0eV whereas the higher 
work function regions have a clear double structure at 
~ 6.0eV and ~ 7.50eV. The structure around 6eV could 
be due to the bulk SiC bands observed along TA, that 
extend from 5.6 to above 6eV^ Other conduction bands 
also disperse along the bulk directions ML and HK. In 
fact, all the high work function spectra show a double 
peak structure in the SE whereas the low work function 
spectra show only one peak. This sort of SE structure 
has already been observed in threshold XPEEM analy- 
sis of graphene on SiC(OOOl). 2 ^ It has been known for 
a long time that bulk graphite produces an intense sec- 
ondary electron signal at 7.5eV above the Fermi level 
(about 3eV above the vacuum),— For example, the pho- 
tocurrent carried by the Bloch constituent of the time- 
reversed LEED wave- function is indeed a maximum near 
7.5eV>22 Thus although the work function difference in 
Fig. [SJa) is small (~ 30meV), we can identify the pres- 
ence of a band at 7.5eV, which for thicker samples, could 
develop into the typical structure of graphite. 



2. Cls core level 

A series of energy filtered images have been acquired 
over the full Cls spectra in the same FoV as the threshold 
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FIG. 7. (a) Intensity map (arbitrary units) of the SiC sub- 
strate component of Cls spectrum after background subtrac- 
tion: the darker regions correspond to lower intensity and 
therefore thicker graphene, FoV53^i. (b) the Cls intensities 
from graphene (squares) and from SiC (diamonds) as a func- 
tion of work function 



data, allowing a pixel-by-pixel extraction of the local Cls 
spectra that can then be directly compared to the work 
function map. The Cls spectrum has two main compo- 
nents, one assigned to the SiC substrate near 283eV, and 
the other to graphene near 285eV. 14 Figure [Ha) shows 
a map of the total area of the SiC component. Because 
thicker graphene regions attenuates the photoelectrons 
from the SiC more, the map in Fig. [TJa) is a good esti- 
mate of the variations in the graphene thickness within 
the FoV. Comparison of the Cls intensity map with the 
work function map in Fig [SVb) shows that there is no 
simple one-to-one correlation between graphene thickness 
and work function. Instead we will show that much of 
the contrast variations are a combination of both work 
function and core level spatial variations, indicating a 
complex chemical structure in the SiC-graphene inter- 
face. 

In order to obtain more chemical-spatial detail, local, 
Cls core level spectra were extracted from each of the 
regions identified in Fig. [7] and are presented in Fig. [5] 
The graphene and SiC intensities show significant vari- 
ations from one region to another, not only in total in- 
tensity but also in terms of the fine structure of both the 
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graphene and SiC components. To fit the spectra, we use 
five peaks, after subtraction of a Shirley background (a 
linear background fit was also tested but does not signifi- 
cantly change the results). Three components are neces- 
sary to fit the main peak of the spectrum that represents 
the graphene. For each of these components a Doniach- 
Sunjic lineshape was used with a 0.2eV Lorentzian width, 
0.3eV Gaussian width and 0.05 asymmetry factor. A one 
component Gaussian lineshape (FWHM of 0.5eV) is used 
to fit the low binding energy peak ascribed to the SiC 
substrate. A Gaussian is more suitable for a wide gap 
semiconductor. The binding energy of this component 
can vary by up to 0.5eV. Finally, a small broad compo- 
nent, which is always present in the Cls spectrum around 
285. 7eV, is ascribed to slight surface contamination due 
to residual gas in the vacuum chamber and fitted using 
a Gaussian with a FHWM = l.leV. The contamination 
component increases slightly for thinner graphene, and 
is largest in region with a work function of 4.34eV. The 
best fits are also reported in Fig. [5] Each spectra is cor- 
related with a region of the sample that is specified by 
its local work function and thickness (determined by it 
Cls as described below). 

The main graphene component lies between 284.30 
and 284.56eV and is always the dominant contribution 
to the spectrum, confirming that there is graphene or a 
graphene like layer over the whole film. This statement 
is supported by the behavior of the SiC component, at 
lower binding energy. It is always much smaller than the 
main graphene component suggesting a near continuous 
graphene coverage. As the graphene intensity increases 
(i.e. the number of graphene layers, increases), the in- 
tensity of the SiC decreases. We can estimate the local 
graphene thickness from the relative SiC and graphene 
Cls intensity. Using a graphene interlayer spacing of 0.34 
and a 1 nm2£ electron mean free path for 654. 3eV pho- 
tons in graphite, the estimated attenuation of the sub- 
strate Cls signal is 22% per graphene layer. Assuming a 
C atom surface densities for graphene (3.8 x 10 15 cm~ 2 ) 
and SiC (1.22 x 10 15 cm~ 2 ), the thickness of the graphene 
within the FoV is between 1 and 3 ML. This estimate is in 
agreement with the spread in the number of ML deduced 
from oscillations in the LEEM backscattered reflectivity 
curves. The thickness values obtained by this method are 
given in Fig. |8l The work function, Cls core level bind- 
ing energies and graphene thicknesses are summarized in 
Table H 

It is immediately obvious that there is not a simple 
one to one correspondence between work function and 
graphene coverage. The same graphene thickness is ob- 
tained in regions with significantly different work func- 
tions. We can group the six different contrast regions in 
Figs. [5jb) and [3(a) into three families. The first consists 
of graphene with 1, 2 or 3 layers. The second, denoted 
by prime superscript, has either 1 or 2 ML graphene, 
while the third, denoted by double prime has only 1 ML 
graphene. 

The complexity of the graphene-SiC interface is re- 
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FIG. 8. Cls core level spectra extracted from the regions 
shown in Fig. [7] together with the best 5-peak fits. The ver- 
tical lines indicate the position of the peak attributed to the 
SiC substrate. The main graphene peak is light grey, the con- 
tamination peak is in black, the SiC component is dark grey. 
The two other graphene peaks, labeled HBE-G and LBE-G, 
flank the main graphene peak. The thickness (calibration is 
described in the text) and work function are given for each 
spectra. 
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TABLE I. Cls BE for the main graphene peak and for the 
SiC peak together, the G/SiC core level intensity ratio for 
different graphene thicknesses and work functions. 



Coverage ^ ( e y) ^ s ( e ^) Ratio 

(ML) WF & ' Graphene Substrate Ig/Is%c 

3 4~45 284.45 282.52 lTX3 

2 4.42 284.55 282.46 10.0 

2' 4.50 284.30 282.84 15.8 

1 4.39 284.56 282.32 5.3 

1' 4.46 284.45 282.75 5.7 

1" 4.34 284.54 282.32 2.0 



vealed by a detailed look at the Cls spectra in Fig. [5] 
There are always two other graphene like peaks besides 
the main peak in all regions, a high binding energy (HBE- 
G) peak around 284. 9eV and a low binding energy (LBE- 
G) peak at 284.0. The intensity of these two peaks is 
lowest for 3 ML graphene and highest for the single layer 
graphene, whatever the family. We therefore associate 
these two structures with carbon below the graphene lay- 
ers, either in the C terminating SiC layer or from some 
of the graphene at the SiC interface. Note that these 
HBE- and LBE-G peak intensities are not the same for 
the three single layer graphene regions. They are most 
intense for the 1" ML region, which also had the strongest 
SE structure associated with the SiC substrate. The 1 
ML and 1" ML graphene have a similar Cls binding ener- 
gies (ABE = lOOmeV). The total work function variation 
over the three single layer graphene regions is 120meV. 
Similarly, 2 ML graphene has significant differences in 
its Cls binding energy and in its work function. The 
main graphene Cls peak BE in the 2' ML graphene is 
shifted 250mcV lower than the 2 ML graphene. At the 
same time, the SiC component shifts 380 meV to higher 
BE in the 2' ML. However, the work function difference 
between the 2 ML and 2' ML regions is 80 meV. The 
work function, Cls graphene, and SiC binding energies 
are plotted together in Fig. O 

It follows that a uniform charge transfer over the whole 
graphene/SiC interface cannot be explained by these ex- 
perimental observations. Charge transfer from the sub- 
strate to the graphene should result in a rigid shift of 
the electronic levels and the work function. However, in 
regions with the same number of graphene layers we ob- 
serve a spread in the work function. The work function 
difference between 1 ML and 1" ML might be ascribed to 
the slightly higher surface contamination, but this can- 
not explain the work function and Cls BE difference with 
the 1 ' ML region which should have the same surface con- 
tamination. The most direct evidence for a non-uniform 
G/SiC interface is the shift of up to 220 meV in Cls 
SiC binding energy for different coverages. The trend in 
the SiC component binding energy emission is shown in 
Fig. If we assume that the SiC Cls signal is dom- 
inated by the first substrate layer, the trend suggests 
two possible interfaces. Indeed, STM studies have shown 




t 1 r 



1 2 3 

Graphene MLs 

FIG. 9. Cls graphene BE (squares), Cls SiC BE (diamonds), 
work function (circles), as a function of graphene thickness. 
The unfilled symbols are for 1, 2, and 3 ML regions, filled 
symbols are for 1' and 2' ML regions. Open symbols with 
dots are for the 1"ML. 

that for very thin UHV grown graphene layers there are 
two different interface reconstructions (2x2 and 3x3)J£ 
The BE variation of the Cls components as a function 
of thickness are also shown in Fig. [9] Different core level 
binding energies for the same graphene thickness further 
supports the interpretation of a non- uniform G/SiC in- 
terface. Comparing the main G and SiC components, we 
see that as the SiC BE increases (charge transfer) there is 
a corresponding decrease in the G BE. Thus, the magni- 
tude of the charge transfer depends on the local interface 
chemistry or structure. 



C. fc-PEEM Results 

In addition to the spatially resolved core level data, we 
are able to create local iso- intensity surfaces in (k Xl k v , E) 
space. In this way, we can immediately visualize the full 
band dispersion in all directions parallel to the graphene 
planes for a selected micron scale region. The objective is 
to correlate the chemical and electronic states obtained in 
XPEEM with a quantitative analysis of the band struc- 
ture of the same micron scale region. Figure [3a) shows 
a complete three dimensional dataset of the band disper- 
sion near the Fermi level of the 2 ML region. The im- 
age series were taken from 2.9eV below the Fermi level 
to 0.3eV above Ep in 50 meV steps, and were repeated 
several times in order to improve statistics without intro- 
ducing camera noise. 

Instead of the single set of six Dirac cones usually re- 
ported for Si-face films^ Fig. HUT a) shows three sets of 
Dirac cones, which we will call principal, secondary and 
tertiary cones. The first two are much more intense than 
the latter. All three have the typical six-fold symmetry, 
although the tertiary cones appear inside the primary 
and secondary cone radius (i.e, at position < TK). Also 
note the symmetry of the primary and secondary cones. 
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FIG. 10. (a) Experimental I(k x ,k y ,E) data collected in in 
the 2 ML region of Fig. [7] The principal and secondary Dirac 
cones at the K and K' points of the first Brillouin zone are 
clearly visible. The secondary cones are rotated by 21.9° with 
respect to the primary cones. Also evident is the already doc- 
umented suppression of the photoemission intensity outside 
the first Brillouin zone due to the pseudo-spin structure when 
using p-polarized light. The tertiary cones are just visible 
inside the double hexagon defined by the principal and sec- 
ondary cones, (b) Schematic of the Brillouin zone of graphene 
showing the two high symmetry directions parallel and per- 
pendicular to TKM' (k x and k y , respectively) extracted from 
the fc-PEEM datasets. An open source volume viewer, de- 
signed for medical imaging, was used to produce the image~i 



They are not full circles because of the suppression of in- 
tensity in the second zone due to matrix element effects 
(an effect well known in graphene) i 32 i 33 However, the ter- 
tiary cones have their symmetry flipped 180°. As we will 
demonstrate below, these three sets of cones are all due 
to two commensurately rotated graphene sheets. 

The secondary Dirac cones have reciprocal lattice vec- 
tors TK rotated by 21.9° with respect to the primary 
cones. This is very close to the value of 21.8° expected 
for a particular pair of commensurate rotated graphene 
sheets^ At first sight, one would be tempted to inter- 
pret the tertiary cones in terms of replicas, like those ob- 
served on the SiC(OOOl) Si-face^ However, those repli- 
cas are due to registry between the (6-\/3x 6-\/3)R30 o re- 
constructed substrate and the graphene overlayer. On 
SiC(OOOl) the graphene layers are known to form com- 
mensurate rotated layers where the average rotation be- 
tween pairs is 30±~7°i£ The tertiary cones in Fig.fTTIarc 
instead due to a diffraction effect caused by the supercell 
formed by the commensurate rotations. To understand 




FIG. 11. Horizontal slice in the 3D fe-PEEM dataset of (a) 2 
ML and (b) 3 ML at a binding energy of 1.3 eV. Secondary 
Dirac cones (graphene reciprocal lattice constant gi) are ro- 
tated by 21.9° with respect to the principal cones (graphene 
reciprocal lattice constant g-z). The tertiary cones are ob- 
tained by diffraction of either the secondary cone by the prin- 
cipal lattice vector g\ (dotted line), or by the primary cone 
by the secondary lattice vector 172 (dotted line). The tertiary 
cones are stronger in 2 ML graphene. 



this, note that the supercell formed by the stacked ro- 
tated sheets is defined by reciprocal lattice vectors Gi 
and Gi- These supercell vectors are, in turn, an integer 
sum of the primary graphene reciprocal vectors, gi a , gib, 
g 2a , and g 2 b, from the two layers (i.e., Gi — pg± a + qg lb 
and Q2 —p'g2a + <l'g2b where p, q, p' and q' are integers). 35 
In Fig. QTJ we show two examples where two of the 
secondary cones are separated by linear sums of lattice 
vectors Gi and G2, i.e., — 3Gi and — (Gi + G 2 )- Note that 
(Gi + G2) — {gia — 32a) indicating that the two rotated 
graphene sheets are indeed a commensurate rotated pair. 
For the purpose of discussion, we use the notation that 
gi 's are the graphene lattice vectors of the primary cones 



10 



k v (A"') k,.(A') M A "> M A '> 




FIG. 12. Band dispersion as a function of k y around the 
7\-point at the first Brillouin zone boundary observed in the 

(a) 2ML and (c) 3ML regions. The respective momentum 
distribution curves over the same energy range are shown in 

(b) and (d). The peak positions are highlighted by the dotted 
lines showing the variation in the position of the Dirac point 
in the different regions. 



and cj2 's are the graphene lattice vectors of the secondary 
cones. Using this notation, it is easy to show that the 
tertiary cones are diffracted replicas of the primary and 
secondary cones. Figure ITT1 shows that a tertiary cone is 
formed by diffracting a primary cone by cj2 or a secondary 
cone by g\ . This diffraction process explains why the ter- 
tiary cones are rotated by 180° relative to the primary 
and secondary cones; they are simply translated Dirac 
cones from a K and K' point of one lattice by the g of the 
second lattice. Because any tertiary cones can be formed 
by diffraction of either of the rotated graphene lattices 
(primary or secondary), the two rotated graphene sheets 
must be stacked on top of each other with a commensu- 
rate rotation. Separated rotated sheets cannot give rise 
to this type of diffraction. 

The 7r-bands of the primary cones cross at a Dirac 
point of 0.30eV for the 2ML film and 0.25eV for the 3ML 
film relative to the Fermi level. Figure [T3] shows cuts 
perpendicular to TKM of the primary and secondary 
cones for 3 ML. The linear dispersion of the 7r-band near 
the Dirac point is evident. Panels (a) and (c) show the 
raw experimental data around the fC-point of the Dirac 
cones. Panels (b) and (d) show the corresponding mo- 
mentum distribution curves (MDCs) extracted from the 
data. The Fermi level calibration was checked by ex- 
tracting the local spectra from a small zone around the 
Dirac cones. This has the advantage of minimizing noise 
and thus allowing a more precise location of the Fermi 
level. The dotted lines show that the 7r-bands cross at 
a Dirac point 50-75meV closer to the Fermi level for the 
principal cones. The group velocity of the quasi-particle 
is 1.0 - 1.1 x 10 6 ms~ 1 in both 2ML and 3ML regions. 

This confirms that the electron doping is a function of 
film thickness and is consistent with the results graphene 
on the Si-face where the charge transfer induced shift of 
the band structure is greater for thinner films^ 



IV. CONCLUSION 

We have studied spatial variations in and correlations 
between the work function, chemical and electronic states 
of few layer graphene grown epitaxially on SiC(OOOl). 
These experiments demonstrate that high thickness uni- 
formity graphene films can be prepared using the CCS 
growth method. In this sample more than 80% of the 
surface is 2-3 layers thick. More importantly, they also 
show that while LEEM can be used to give local infor- 
mation on the graphene films thickness, LEEM contrast 
variation on C-face graphene films are not simply due 
to film thickness alone. X-PEEM's spatial resolution has 
proven to be very useful in illuminating the origin of these 
contrast variations. The O.l^m spatial resolution in the 
present experiment is much better than that currently 
obtained in area averaged high resolution ARPES with 
synchrotron radiation. The work function variations in 
Fig. derived from the X-PEEM, combined with the lo- 
cal Cls spectra in Fig [5] show an interface chemistry that 
is more complex than originally suspected. 

To demonstrate this, Fig. [13] shows an area averaged 
Cls spectrum from the 53/im FoV. The overall spec- 
trum is broader and the fine structure visible in the local 
spectra of Fig. [3] is smeared out in the area averaged 
spectrum. Without the spatial resolution, one would 
have concluded that only one broad substrate peak was 
present rather than a narrower peak whose energy posi- 
tion shifts depending on the graphene thickness. Simi- 
larly, the area averaged work function masked the com- 
plexity of charge transfer at the graphene-SiC interface. 
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FIG. 13. Area averaged Cls core level spectrum. The average 
is over the 53/^m FoV 

While the X-PEEM data agrees with the general trend 
observed by others that the work function trends to 
higher values with thicker graphene films, it shows that 
the situation is more complicated. For a given thickness, 
both the local work function and Cls BE vary apprecia- 
ble. The data suggests that at least two SiC terminations 
with different local bonding appear to be present at the 
interface. We note that the long range order of these 
different bonding areas is not very high since the micro 
LEED never sees any superlattice diffraction spots. For 
a given graphene-SiC interface, the core level data can 
be interpreted within the framework of a charge trans- 
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fer model between the substrate and the graphene that 
rigidly shifts the Cls core level. The effect of the charge 
transfer on binding energy is greater for carbon closer to 
the interface, i.e. for carbon in thinner graphene films. 
The correlation between the work function and the core 
level BE strongly supports this interpretation. On the 
basis of these results we suggest that there are two dis- 
tinct interfacial chemical regions. 

It is possible that these two regions are similar to 
the 2x2 and 3x3 observed by STM during the early 
UHV growth of C-face graphene^ In the furnace grown 
graphene used in this results, the higher growth tem- 
perature (1550C versus 1100C) prevents any significant 
interface ordering so that no LEED spots are observed. 
However, Hiebel et ali£ suggested that there is a stronger, 
localized interaction between the 2x2 reconstructed C 
layer and the SiC. Thus, one would expect a higher den- 
sity of sp 3 C-Si bonding between the C termination layer 
and the underlying SiC. From simple charge transfer ar- 
guments, this would translate into a lower Cls binding 
energy. On this basis the Cls spectra of the V-2' ML 
regions could be tentatively associated with a local 2x2 
surface. In that case regions labeled as 1, 2 and 3 ML in 
Table Q] would be associated with a local 3x3 interface. 
However, it is also possible that these different regions 
of the surface may be due to local silicon concentrations. 
This can only be resolved by performing Si 2p mapping 
(experiments that are now being planned). 



We have also demonstrated that fc-PEEM gives de- 
tailed band structure information from commensurately 
rotated graphene sheets. The easy switching between real 
and fc-space imaging modes allows, within the limits of 
the spatial resolution and field aperture size, correlation 
between the chemical and electronic states of the surface 
and graphene-SiC interface. In this study we conclusively 
show that commensurately rotated graphene pairs exist. 
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